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Abstract 



We explore generic ground-state and low-energy statistical properties of 
many-body bosonic and fermionic one- and two-body random ensembles 
(TBRE) in the dense limit, and contrast them with Random Matrix Theory 
(RMT). Weak differences in distribution tails can be attributed to the regu- 
larity or chaoticity of the corresponding Hamiltonians rather than the particle 
statistics. We finally show the universality of the distribution of the angular 
momentum gap between the lowest energy levels in consecutive J-sectors for 
the four models considered. 
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Wigner introduced the Gaussian orthogonal ensemble (GOE) to deal with the statistics 
of high- lying levels of many-body quantum systems . Although the fluctuations of certain 
observables predicted by the GOE agree well with experimental observations, it represents a 
system in which all the particles interact simultaneously and this is a priori not appropriate 
to describe many-body systems in which the two-body interaction is predominant 
The two-body random ensemble (TBRE) was introduced to improve upon the physical 
limitations of RMT [Q. The TBRE Hamiltonian includes only up to two-body operators, 
whose coefficients are real random numbers. This ensemble reproduces a Gaussian level 
density and the GOE level repulsion as desired, but until recently, only the dilute limit was 
analytically tractable 0. This limit corresponds to having a large number of particles {Np 3> 
1) and an even larger number of single-particle states K ^ 1, so that Np/K <^ 1. In such 
a limit Pauli's principle has only a marginal effect, so particles statistics are unimportant. 
Many physical regimes lie outside the dilute limit: for bosonic particles, we may consider a 
large number of particles with a finite K -C A''^; similarily, fermionic antisymmetry becomes 
relevant as one approaches half-filling, Np ^ K/2. Since we are interested in those aspects, 
we do not restrict ourselves to the dilute limit. It has only recently been shown, from the 
connection between the Lanczos tridiagonalization and random polynomials, that this case 
can also be treated analytically in some instances |0. 

We investigate the statistical properties of the low-lying levels of interacting complex 
systems, extracting from the edge of the spectra, generic properties that would refiect the 
basic structure of many-body systems such as molecules, atomic nuclei or quantum dots close 
to, or in their ground-state. We compute the lowest eigenvalue distribution for two bosonic 
and two fermionic TBREs and show that, after a proper rescaling, the distributions have 
a surprisingly weak dependence on the nature of the particles. Comparing to predictions 
from RMT lead to surprising similarities. The fact that the distribution of the lowest energy 
state depends only weakly on the nature of the particles and not on the ensemble used to 
study it points to its universality, but a closer inspection of the tails of these distributions 
indicate that they depend on the chaoticity of the system analyzed. To understand the role 
of chaos, we study models which are regular as well as chaotic and find that the tails of the 
ground state distributions for chaotic systems are closer to that obtained from the large-N 
GOE, while for integrable systems it is closer to the N=2 GOE. The distribution of energy 
differences between the lowest energy levels in two consecutive angular momentum sector 
(J and J + 1) is found to be more robust, having the same shape in all four cases within 
numerical accuracy (deviations occur for the single j-shell model at small gap values). Since 
this distribution depends on correlations between Hilbert subspaces with different quantum 
numbers, it is outside the scope of RMT. We expect that this distribution is generic, a 
conclusion which is borne out by the bosonic nature of the corresponding excitation. Our 
analysis is based on numerical and analytical results. The analytical treatment is now 
possible not only for the TBRE 0, as already pointed out, but also for the large N limit of 
the GOE, for which an expression for the distribution function of the largest eigenvalue was 
obtained in terms of a particular Painleve II function 

We start by summarizing analytical results obtained for the GOE. A general matrix H of 
this ensemble can be expressed in terms of a rotation matrix O and the diagonal eigenvalue 
matrix E = diag(i?i, E2, ■■■Em) as H = O^EO. The distribution of the ground state energy 
is P{Eq) = {6{Eq — El)), where (■ ■ ■) indicates GOE averaging and the probability of having 
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a matrix element Hij is P{Hij) cx exp(— |Trif^). Two limits are considered: N = 2 and 
oo. For N = 2, direct integration gives 



^2^°^(^o) = ^|exp(-Eo^)-y|exp(-EoV2)i?o erfc 



Eq 



(1) 



This distribution is not readily calculated for arbitrary N, however in the limit N oo, 
Tracy and Widom derived an expression for the distribution of the largest eigenvalue ( 
1^), which in terms of the lowest eigenvalue is written as 



P^^^(Eo) = - [^J_^^ q{xfdx + g(-^o) j exp + EM^fdx + q{x)dx 



(2) 



Here g(-Eo) satisfies the Painleve II equation and the boundary condition is g(-E'o) ~ Ai(_E'o) 
a.s Eq ^ oo, where Ai is the Airy function. The numerical solutions for this expression are 
plotted on the top panel of Fig.l, where one sees that it agrees well with the distribution 
for a GOE of A^ = 200. This function has the asymptotics 

\ogPg^'^{E,)^-\Eo\\ Eo^-oo, logP^«^(Eo)~-i?o'^ Eo oo. 

(3) 

In contrast to A^ = 2, the tails deviate from Gaussian behavior. For completeness, we also 
show the average ground state energy for the GOE, although this quantity is not expected to 
have physical consequences. We find that (£"0) = 2.01{2)y/N + 1.58{2)N-^/^, which agrees 
with the predicted functional form ^j^. On the bottom panel of Fig.l, we show that this 
expression fits numerical results for A^ G [2, 1000]. 

In contrast to the GOE, the TBRE is not a purely statistical model, but rather is 
constructed in a particular bosonic or fermionic space. We take it in the form 

H = Y1 (^ijO'tO'j + Vijkiafajakai. (4) 

ij ijkl 

When originally introduced, the Hamiltonian operators of the TBRE had only (interaction) 
two-body terms, but we include the one-body term, a^(aj) represents boson or fermion 
creation (annihilation) operator and the coefficients eij and Vijki are taken as Gaussian 
random variables once certain physical constraints are imposed such as rotation invariance, 
time reversal invariance and conservation of total spin, isospin and parity. 

As already mentioned, some analytical results for bosonic |^,^ and fermionic 0,^ models 
in the dilute limit were obtained, but here we consider the non-dilute limit, for which little 
is known. We want to find out how this limit and some other model dependencies, such as 
the nature of the particles or the system's integrability, may affect the statistical properties 
of many-body systems. 

Bosonic Models 

We study two bosonic models, the U(4) vibron and U(6) interacting boson models, used 
to explore molecular or nuclear collective excitations. In the U(4) TBRE model, bosons 



with J"^ = O"*" [s, s~^) and J'^ = 1 {p, p^) are coupled to form a scalar Hamiltonian |jTO| 
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Hu(4) = (esS^s + epP^.p)/Np (5) 
+ {j[p^P^f^-[ppf^ + j[P^P^f^-[Ppf^ + ^[s+s+f\[ssf^ 
+ y + ^{[s^s^f\[ppr + h.c.)}/Np{Np - 1). 

The square brackets denote angular momentum coupling, the dots represent scalar products 
and Np is the total number of bosons. Since the 1- and 2-body matrix elements are propor- 
tional to Np and Np{Np — 1) respectively, scaling allows all coefficients in (|]) to be Gaussian 
random numbers of unit variance. For the choice of coefficients above, the lowest eigenvalue 
distribution has been obtained analytically, and in the large Np limit 

P{Eo) oc e"°-^^oerfc(0.663Eo) + 0.423e-^o/2-^^erfc(0.632£'o). (6) 

The second model we investigate is the nuclear U(6) model (IBM) which consists of scalar 
J'^ = 0"*" (s, s"*") and quadrupole J'^ = 2+ {d, d'^) bosons coupled by one- and two-body 
interactions. The Hamiltonian has the form |TI[] 

Hui6) = (e^s+s + tdd'^.d)/Np 

L=0,2,4 

+ a^[s^d^]^^\[sdf^ + a4{s^ s^f\[ddp + h.c.}/Np{Np - 1), (7) 

We will consider Np = 16, which is typical for a heavy collective nucleus, and the distribution 
for the ground state energy is obtained numerically. An important distinction between the 
U{4) and U{6) models is that the former is integrable for all parameters, while the latter 



is generally chaotic ||I2[. Chaos usually applies to fluctuation properties near the middle of 
the spectrum, rather than at the edges. Never the less, there is a subtle difference between 
these models which we attribute to the underlying chaos, as we discuss below. The ground 
states of these models also have well known phase transition behaviors. In spite of this, the 
ground state distributions will be seen to be surprisingly insensitive to this. 
Fermionic Models 

The two fermionic models we investigate are a model for interacting nucleons in a single 
j-shell and a model for randomly interacting electrons. For the single j-shell we use Np = 6 



particles in the j = 15/2 shell. The matrix elements of the Hamiltonian are |]T3 



< f^aJM\ ^ Vik\3^^a'JM > 

i<k 

^Np{N^ E [j^^«J{|j^-^(«V'Oj^(J0^][j^^-'(«V'')j^(J')^|}j^''«'J]t;,,, (8) 

^ a", J", J' 

where J is the total angular momentum, M is the z-projection, a denotes all other quantum 
numbers, the term in brackets corresponds to the coefficients of fractional parentage and vji 
are random numbers. 

The second fermionic model is a model for randomly interacting j = 1/2 fermions (lEM) 
and is given by 
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where the one-body spectrum is Wigner-Dyson distributed with G [—K/2;K/2] {K/2 
gives the number of spin-degenerate orbitals, the average level spacing is then A = 1), 
U2p G [—U; U] are Gaussian random numbers and j(j') =T (i) are spin indices. We study 
up to A'p = 7 fermions on K/2 = 7 orbitals. 



While the single j-shell model describes non-chaotic system , the spectral properties 
of the lEM depend on the ratio of the interaction strength to the one-body level spacing 
U/A: for f//A -C 1/{KN^) the model is non-chaotic with a Poissonian spectrum while for 



U/A^ l/(KNp) it is chaotic with a Wigner-Dyson distributed spectrum |jT6|. 

To compare the lowest eigenvalue distribution for these four models with the ones for the 
GOE, we re-center the distributions and rescale their widths by using the following scaling 
prescription 

eo = . (10) 

This removes model dependencies associated with scales, level densities and so forth. We 
show in Fig. 2 distributions of the rescaled ground state energy eo for the four models de- 
scribed above. The figure splits into regular systems on the top two panels and chaotic 
ones on the two bottom panels. Bosonic models are left and fermionic ones are right. The 
chaoticity of the lEM depends on U/A and the data shown on the lower right panel of Fig. 2 
for the lEM correspond to ?7/A = 1, well into the chaotic regime |jl6| . 



One clearly sees that the ground state energies have approximately the same distribution 
for the four TBREs and the GOE of small and large dimension - this seems to be a robust 
property, and in particular, these distributions seem not to depend on the nature of the 
particles. Moreover, the agreement of the distributions with the ones obtained from the GOE 
demonstrates that the distribution of the ground state energy is also not sensitive to the type 
of ensemble used. A closer inspection of the presented data does reveal discrepancies between 
different models in the distribution tails. The lowest energy distribution for regular systems 
coincides better with the distribution for the two-dimensional GOE, while the distribution 
for chaotic systems is in better agreement with the one obtained with a large-N GOE. The 
GOE dimension necessary to describe a system can be associated with the integrability of the 
model. This conclusion is corroborated by the data shown on Fig. 3 for the lEM at various 
U/A, both in the integrable and chaotic regimes. Evidently the crossover from Poisson to 
Wigner-Dyson distributed level spacing (i.e. from integrable to chaotic quantum dynamics) 
is accompanied by a crossover from two-dimensional to infinite-dimensional GOE fitting 
of the ground-state distribution. We also found that the ground-state distribution for the 
lEM exactly at half-filling {K = 14 and Np = 7), gets close to a two-dimensional GOE, as 
shown at the top of Fig.4. Since the associated spectrum is still Wigner-Dyson distributed, 
we conclude that this is a manifestation of the particle statistics. The only influence of 
the latter therefore is to single out half-filled fermionic systems which despite their chaotic 
character have a distribution of ground-state energies which is closer to the N = 2 GOE. 

Consider now distributions which cannot be addressed through RMT, namely distribu- 
tions which involve correlations between Hilbert subspaces with different quantum numbers. 
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One quantity of fundamental interest is the energy gap. In the models we study, this typi- 
cally involves states of different spin: for the vibron model it is e{J'^ = 1~) — e(J'^ = O"*"), 
for the IBM 6{.r = 2+) - e{.r = 0+), while for the single j-shell model and the lEM 
e{J = 1) — e{J = 0). We will call this energy difference the spin gap. The ground state 
is mostly dominated by 0^ (this predominance almost reaches 100% for the lEM [|I^) and 
the lowest level with a different spin usually corresponding to the first excited state (this is 
however not the case for the lEM), so we are looking at the very edge of the spectrum. Even 
though we have both bosonic and fermionic models, the distributions are very similar and 
also agree very well with the analytical expression derived for the vibron model ^ as it is 
evident from Fig. 5. The distribution of the spin gaps is thus also a generic property which 
we attribute to the bosonic nature (i.e. 6 J is an integer) of the corresponding excitation in 
all four models studied here. This explains why the particle statistics have no influence here. 
Note that in the case of the lEM, universality of the spin gap requires a finite interaction 
strength U/A ^ 1. In the limit U/A 0, the spin gap is determined by the one-body 
spectrum ea and thus given by a GOE one- (two-) spacing distribution for Np even (odd). 

In summary, investigating various TBREs, we have found that in spite of their funda- 
mental differences, they have strikingly similar features in both their ground state energy 
distribution and their spin gap distribution. The ground state energy distributions have 
common features, and coincide to GOE distributions of either small or large dimensions, 
for regular or chaotic TBREs respectively. We conclude that the edge of the spectra ex- 
hibit two unexpected generic properties: the distribution of the ground state energy and the 
distribution of the energy differences between the lowest levels with different spins. This 
latter point is presumably related to the apparent regular structure found in the angular 
momentum structure of low-lying levels for several TBREs [|l^] . 
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FIGURES 




FIG. 1. Top: Analytical distribution of the ground state energy for GOE of dimension N = 2 
(solid line) and N = oo (dashed line) . The numerically obtained distribution for N = 200 (squares) 
already coincides with the N = oo GOE. Bottom: Average GOE ground-state energy {Eq) as a 
function of the matrix size N. The solid line gives the theoretical estimate. 
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FIG. 2. Distributions of the ground state energy for the vibron model (top left), the single 
j-shell model (top right), the IBM (bottom left) and the lEM (bottom right) for Np = 3 (circles) 
and 4 (squares). In all panels, solid and dashed lines give the distribution for the N = 2 and 
N = oo GOE respectively. The top panels correspond to non-chaotic systems, for which the data 
are well fitted by a = 2 GOE. The bottom panels correspond to chaotic systems, where data 
are closer to the N = oo GOE. 
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FIG. 3. Crossover to quantum chaoticity for the lEM with Np = A and U/A = 10~^ (diamonds), 
0.05 (triangles), 1 (squares) and oo (circles). Symbols for data in the integrable regime are filled, 
while solid and dashed lines give the distribution for the N = 2 and N = oo GOE respectively. 
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FIG. 4. Top panel: Distribution of the ground state energy for the lEM for Np = 4 (triangles), 
5 (crosses) and 6 (stars) . Bottom panel: Distribution of the ground state energy for the lEM at half 
filling, Np = 7 (circles). In both panels, solid and dashed lines give the distribution for the N = 2 
and N = oo GOE respectively. The particle statistics plays a role at half-filling only, rendering the 
distribution closer to the N = 2 GOE. 
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FIG. 5. Spin gap distribution for the U(4) vibron model (x) the U(6) IBM (+) the single 
j-shell model (•) and the lEM with Np = 5 (triangles), 6 (squares) and 7 (diamonds). The solid 
line corresponds to the analytical expression obtained for the vibron model. The only significant 
deviations occur for the single j -shell model for low gap values. 
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